Microscopic theory of spin-triplet /-wave pairing in quasi- ID organic superconductors 
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We present a microscopic theory of fluctuation-mediated pairing mechanism in organic supercon- 
ductors (TMTSF) 2 X, where the experimentally observed coexistence of 2kp charge fluctuation and 
2k,F spin fluctuation is naturally taken into account. We have studied, within the random phase 
approximation, the extended Hubbard model at quarter-filling on a quasi-one-dimensional lattice, 
where we consider the off-site repulsive interaction up to third next nearest neighbors along with 
the on-site repulsion. The results show that spin-triplet /-wave-like pairing can be realized in this 
system, dominating over singlet ci-wave-like pairing, if 2fcp spin and 2kF charge fluctuations coexist. 



It has been a long standing issue to clarify the super- 
conducting state of quasi-one-dimensional (qlD) organic 
superconductors (TMTSF) 2 X {X = PF 6 , C10 4 , etc.), so 
called the Bechgaard salts^i. One may expect uncon- 
ventional superconducting states due to the quasi-one- 
dimensional nature of these materials as well as the elec- 
tron correlation effects generally seen in organic materi- 
als. In fact, there have been experiments suggesting the 
existence of node a 3 i 4 i 5 , while recent experiments showing 
a large H c ^^ and an unchanged Knight shift across T<A£ 
supports a realization of spin-triplet pairing. The pairing 
symmetry of (TMTSF) 2 X has thus become quite a hot 
issue very recently. 

Theoretically, triplet p-wave pairing state in which the 
nodes of the pair potential do not intersect the Fermi sur- 
face has been proposed from the early days&ifiiAI. How- 
ever, the occurrence of triplet pairing in (TMTSF) 2 X 
is curious from a microscopic point of view since su- 
perconductivity lies right next to a 2kp spin den- 
sity wave (SDW) in the pressure-temperature phase 
diagram 12 . Naively, SDW spin fluctuations should favor 
spin-singlet d-wave-like pairing as suggested by several 
other authorsi^iii^. However, one should note that the 
insulating phase is not pure SDW at least for some an- 
ions, namely, 2kF charge density wave (CDW) actually 
coexists with 2kp SDW In fact, one of the present 
authors has phenomelogically proposed 2 ^ that triplet f- 
wave-like (see Fig. 1(c) for a typical pair potential) pair- 
ing may dominate over p-wave pairing and become com- 
petitive against d- wave-like pairing (see Fig. 1(b)) due 
to a combination of quasi-lD (disconnected) Fermi sur- 
face and the coexistence of 2fcp SDW and 2fcp CDW 
fluctuations. A similar scenario has been proposed by 
Fuseya et alw* Concerning the /-wave vs. d-wave com- 
petition, it has also been proposed that magneto tun- 
neling spectroscopy 2 ^ via Andreev resonant states^ is a 
promising method to detect the /-wave pairing. 

However, there has been no microscopic theory for 
/-wave pairing in (TMTSF^X starting from a Hamil- 
tonian that assumes only purely electronic repulsive 
interactions^* To resolve this issue, here we study an ex- 
tended Hubbard model at quarter-filling on a quasi-one- 
dimensional lattice within the random phase approxima- 
tion (RPA). We consider off-site repulsions up to third 



nearest neighbors along with the on-site repulsion in or- 
der to naturally take into account the coexistence of 2kp 
charge and 2kp spin fluctuations. The merit of adopting 
RPA 25 is that we can easily take into account the off-site 
repulsion as compared to fluctuation exchange (FLEX) 
approximation, where it is by no means easy to take into 
account distant interactions 2 ^^. 
The model Hamiltonian is given as 



H 



E 

<l,j>,<7 



tij C icr C jcr' 



» \i—j\= m l 



where c irT creates a hole (note that (TMTSF^X is actu- 
ally a 3/4 filling system in the electron picture) with spin 
a =t, | at site i=(i a ,if,). Here, < i,j > stands for the 
summation over nearest neighbor pairs of sites. As for the 
hopping parameters, we take £y = t a for nearest neigh- 
bor in the (most conductive) a-direction, and = tf, for 
nearest neighbor in the 6-direction. We choose % = 0.2£ a 
to take into account the quasi-onc-dimcnsionality. t a is 
taken as the unit of energy throughout the study. U 
and Vij are the on-site and the off-site repulsive inter- 
actions, respectively. We take distant off-site repulsions 
because it has been shown previously that nearest neigh- 
bor and second nearest neighbor off-site repulsion is nec- 
essary to have coexistence of 2k p spin and 2k p charge 
density wavesiiSiia Here we take off-site repulsions up to 
third nearest neighbors, namely, Vi,j — Vq, V\ and V% 
with m=\ i a — j a | = 1, 2 and 3, respectively. The effect of 
the third nearest neighbor repulsion, V2 , will be discussed 
at the end of the paper. The effective pairing interactions 
for the singlet and triplet channels due to spin and charge 
fluctuations are given as 

V s (q,a Jl )^U + V(q) + ^U 2 Xs (q,u; l ) 



--(U + 2V(q)) 2 Xc (q,LU l ) 



V t (q,LU l ) = V(q)--U 2 Xs (q,LU l ) 



-(U + 2V(q)) 2 X c(q,0J l ) 



(1) 



(2) 
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within RPA, where 

V(q) = 2V cos q x + 2Vi cos(2fe) + 2V 2 cos(3g x ) (3) 

and loi is the Matsubara frequency. Here, Xs and Xc are 
the spin and charge susceptibilities, respectively, which 
are given as 



xo(q,ui) 

1 - Uxo(q,ui) 



l+(U + 2V(q)) X a(q,ui) 
Here xo is the bare susceptibility given by 



(4) 



(5) 



Xo(<?,w z ) 



1 E 



/( e P+q) - f( e p) 



with £fe = — 2t a cosk a — 2tbCoskb — [i and f(e p ) = 
l/(exp(e p /T) + 1). xo peaks at the nesting vector Q 
(= (7r/2,7r) here) of the Fermi surface. The terms pro- 
portional to Xs and Xc in eqs. Q and @ represent 
effective pairing interactions due to the spin and charge 
fluctuations, respectively. The chemical potential fi is 
determined so that the band is quarter-filled, which is 
/i = —1.38 here. We take U = 1.6 throughout the study, 
which is large enough to have strong 2kp spin fluctua- 
tions (large Xs{Q) = Xs{Q,0)) but not so large as to 
drive SDW instability at high temperatures. The off-site 
interactions are chosen so that 2fcp charge fluctuations 
are induced as discussed later. In the actual numerical 
calculation, we take N = 400 x 40 fc-point meshes ex- 
cept for low temperatures, where we take N = 800 x 80 
meshes. 

To obtain the onset of the superconducting state, we 
solve the the gap equation within the weak-coupling the- 
ory, 



AA(fc) = -^F^ t (fc-fc',0) 



tanh(/3e k //2) 
2^ 



A(fc'). (6) 



The transition temperature Tc is determined by the con- 
dition, A = 1. In the weak coupling theory, lo dependence 
of the pair potential A(fc) is neglected. Although this ap- 
proximation is quantitatively insufficient, it is expected 
to be valid for studying the pairing symmetry of A(fc) 
mediated by both spin and charge fluctuations. In the 
following calculations, we study triplet and singlet cases 
with A(fe) = -A(-fc) and A(fe) = A(-fe), respectively. 
We define <j) s (k) = A(fe)/A M and <f> t (k) = A(k)/A M for 
singlet and triplet pairing, respectively, where Am is the 
maximum value of the pair potential. 

From eqs.J*[J-((SJ|, it can be seen that when U ~ — {U + 
2V{Q)) is satisfied, \V S {Q)\ ~ |V*(Q)| holds, apart from 
the first order terms, which are negligible in the limit of 
strong spin and/or charge fluctuations (but turn out to 
be important in the actual cases considered later). This, 



along with the disconnectivity of the Fermi surface (note 
that the number of nodes intersecting the Fermi surface 
is the same between /- and d- waves) , is expected to make 
spin triplet /-wave pairing competitive against singlet d- 
wave pairing^* 

We now move on to the calculation results. First, we 
focus on the case where spin fluctuation is dominant, e.g., 
Vq = V\ = Vi = 0. As shown in Fig. 1, the magnitude 
of A for the singlet case is much larger than that for the 
triplet case. The resulting singlet pair potential 4> s (k) 

changes sign as H 1 — along the Fermi surface (see Fig. 

1(b)). We call this d-wave pairing, where s (fe) is roughly 
proportional to cos(2fe x ). On the other hand, the triplet 

pair potential 4>t(k) changes sign as H 1 1 — along the 

Fermi surface (see Fig. 1(c)). We call this /-wave, where 
4>t(k) is roughly proportional to sin4fc a . The results here 
are expected from the previous FLEX studyfSS, 

Let us now move on to the cases where we turn on 
the off-site repulsions. In order to to have the coexis- 
tence of 2kp spin and 2kp charge fluctuations as exper- 
imentally observed, namely, to have Xs(Q) ~ ~Xe(Q) = 
— Xc(Q)0), U ~ — (U+2V(Q)) has to be satisfied as men- 
tioned earlier. To accomplish this, V\ has to be close to 
Z7/2, as can be seen from eq.(|3J|. Namely, since the x com- 
ponent of Q is Q x = it/2, the V\ term in eq.@ is domi- 
nant for q~Q, making U - -(U+2V(Q)) if V 1 ~ U/2. 
Thus we first choose Vi — U/2 — 0.8. As for the other 
off-site repulsions, we first choose Vb = 1.2 and V<x = 0.5 
as a typical value. As shown in Fig. 2(a) and 2(b), 
Xs and Xc peak around {k a ,k b ) = ±(n/2, n),±(n/2, -n) 
and the maximum values are both about 9.1, so that we 
indeed have the situation where 2fcp spin and 2kp charge 
fluctuations coexist. As shown in Fig. 2(c), the magni- 
tude of A for triplet pairing is now much larger than that 
for singlet pairing. The corresponding singlet pair po- 
tential 4> s {k) has the d-wave form as shown in Fig. 2(d), 
while the triplet one (f>t(k) has the /-wave form as shown 
in Fig. 2(e). Note that the result of Applet A s i ng i e t 
is rather unexpected from the previous phenomelogical 
arguments . because when Xc(Q) = Xs(Q), /-wave is only 
degenerate with d-wave in the previous theory. The origin 
of this discrepancy is the first order terms in eqs.Q and 
(J2J, which are neglected in the phenomelogical theory. 
Thus, we have obtained a remarkable result here, namely, 
/-wave can completely dominate over c?-wave when 2kp 
spin and charge fluctuations coexist. 

In order to further look into this point, we next reduce 
V\ from 0.8, thereby suppressing the charge fluctuation. 
The maximum value of Xc (not shown) is 5.2 and 3.2 for 
Vi = 0.78 and Vi = 0.75, respectively. Although the 
maximum value of Xc is smaller than that of Xs (= 9.1) 
in these cases, A for the triplet case is still larger than 
(Vi = 0.78) or competitive against {V\ = 0.75) that for 
the singlet case as seen in Fig. 3 . This means that /-wave 
pairing has a chance to be realized even if spin fluctuation 
dominates, as far as 2kp charge fluctuation exist to some 
extent^ 

Finally, in order to look into the effect of the third 
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nearest neighbor interaction V2, we set V2 = leaving 
the other parameters the same as in Fig. 2. 

As seen in Fig. 4(a), the peak of Xc is slightly shifted 
towards Akp (q x = tt) when V2 = 0, compared to the 
case of V2 = Q-5t shown in Fig. 2(a). (xs is the same as 
Fig. 2 (a) since V2 does not affect Xs within the present 
formalism.) In this case, the singlet-triplet competition 
becomes much more subtle as seen in Fig. 4(b). The rea- 
son for this can be found in Fig. 4(c), namely, <p s (k) in this 
case has the same sign on most portion of the Fermi sur- 
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FIG. 1: Calculation results for U = 1.6, Vb = Vi = V2 = 0. 
(a) Temperature dependence of A for singlet (solid line) and 
triplet (dotted line) pairings. Contour plots of (b) cj> s (k) and 
(c) 4>t{k) at T = 0.01. In (b), and (c), the solid lines represent 
the Fermi surface and the dotted lines denote the nodal lines 
of the pair potentials. 
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FIG. 2: Calculation results for U = 1.6, Vb = 1.2, Vi = 0.8, 
Vs = 0.5: (a) X s(k) at T = 0.01. (b) Xc (Jfe) at T = 0.01. 
(c) Temperature dependence of A for singlet (solid line) and 
triplet (dotted line) pairings. Contour plots of (d) 4> s {k) and 
(e) 4>t(k) at T = 0.01. In (d), and (e), the solid lines represent 
the Fermi surface and the dotted lines denote the nodes of the 
pair potentials. 
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FIG. 3: Calculation results for U = 1.6, Vo = 1.2, and 
V2 = 0.5. (a)Temperature dependence of A for singlet 
(solid line) and triplet (dotted line) pairings for V\ = 0.78. 
(b)Temperature dependence of A for singlet (solid line) and 
triplet (dotted line) pairings for V\ = 0.75. 



face. In other words, it is more like s-wave than d-wave. 
This 's-wave' pairing is induced by charge fluctuation, 
which does not totally cancel out with spin fluctuation in 
eq. 0J) because the wave vector at which Xc peaks deviates 
from that for Xs- Since the 's-wave' pair potential has the 
same sign on most portion of the Fermi surface, almost all 
the pair scattering processes on the Fermi surface, medi- 
ated by the attractive interaction (note the minus sign in 
eq.Q) due to charge fluctuation, have positive contribu- 
tion to superconductivity, making singlet pairing much 
more enhanced compared to the case with nonzero V2- 
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FIG. 4: Calculation results for U = 1.6, Vo = 1.2, Vi = 0.8, 
V2 = 0: (a) Xc(fe) at T = 0.01. (b)Temperature dependence 
of A for singlet (solid line) and triplet (dotted line) pairings. 
Contour plots of (c) s (fc) and (d) cj> t (k) at T — 0.01. In (c) 
and (d), the solid lines represent the Fermi surface and the 
dotted lines denote the nodes of the pair potentials. 



Conversely, the present results show that V2 has the ef- 
fect of stabilizing 2kp charge fluctuation, which has a 
tendency to shift towards 4/cf fluctuation when only Vo 
and V\ are present, and this effect in turn suppresses sin- 
glet pairing because in that case a strong cancellation 
occurs between the third and the fourth terms in eq.Q. 
Since the screening effect is known to be weak in quasi- 
one-dimensional systems, it is likely that such a distant 
off-site repulsion is present in the actual (TMTSF^X. 



To summarize, we have presented a microscopic the- 
ory of pairing mechanism in organic superconductors 
(TMTSF) 2 X , where we have taken into account the coex- 
istence of the 2kp charge fluctuation and 2kp spin fluctu- 
ation by considering off-site repulsions up to third nearest 
neighbors. We have shown that the /-wave triplet pair- 
ing symmetry can be realized in this system when 2kp 
charge density fluctuation and 2kp spin density fluctu- 
ation coexists. Surprisingly, the condition for realizing 
/-wave pairing is much more eased compared to that in 
the previous phenomelogical theorySl. 



In the present paper, we have neglected the realis- 
tic shape of the Fermi surface observed in the actual 
(TMTSF) 2 X2£. Although the influence of this effect on 
the /-wave pairing is expected to be small because the 
x-component of the nesting vector is close to tt/2 in any 
case, a detailed analysis remains as a future study. Ex- 
perimentally, it would be interesting to verify the real- 
ization of /-wave superconducting state by some phase 
sensitive probesiS&Si 
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